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ABSTRACT 

Planetary companions to the source stars of a caustic-crossing binary mi- 
crolensing events can be detected via the deviation from the parent light curves 
created when the caustic magnifies the star light reflecting off the atmosphere or 
surface of the planets. The magnitude of the deviation is 6p ~ epPp^^^, where ep 
is the fraction of starlight reflected by the planet and pp is the angular radius 
of the planet in units of angular Einstein ring radius. Due to the extraordi- 
narily high resolution achieved during the caustic crossing, the detailed shapes 
of these perturbations are sensitive to fine structures on and around the plan- 
ets. We consider the signatures of rings, satellites, and atmospheric features on 
caustic-crossing microlensing light curves. We find that, for reasonable assump- 
tions, rings produce deviations of order 10%5p, whereas satellites, spots, and 
zonal bands produce deviations of order l%6p. We consider the detectability of 
these features using current and future telescopes, and find that, with very large 
apertures (>30m), ring systems may be detectable, whereas spots, satellites, and 
zonal bands will generally be difficult to detect. We also present a short discus- 
sion of the stability of rings around close-in planets, noting that rings are likely 
to be lost to Poynting-Robertson drag on a timescale of order 10^ years, unless 
they are composed of large cm) particles, or are stabilized by satellites. 
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Subject headings: gravitational lensing - planets: rings - planets and satellites: 
general 



1. Introduction 

Precise radial velocity surveys have detected over 100 planetary companions to FGKM 
dwarf stars in the solar neighborhood (see http://cfa-www.harvard.cdTi/planets/catalog for 
a list of planets and discovery references). Among the interesting trends that have been 
uncovered in this sample of planets arc a positive correlation between the frequency of planets 
and metallicity of the host stars (Gonzalez 1997, 1998; Laughlin 2000; Santos, Israelian, & 
Mayor 2001; Reid 2002), a paucity of massive, close-in planets (Zucker & Mazeh 2002; Patzold 
& Rauer 2002), and a 'pihng-up' of less- massive, close-in planets near periods of P ~ 3 days. 
This latter trend is important because the number of planets which transit their parent 
stars is roughly proportional to 1/a, where a is the semi- major axis. The discovery and 
interpretations of these global trends provide clues to the physical mechanisms that affect 
planetary formation, migration, and survival. 

A somewhat different way of obtaining clues about the physical processes at work in 
planetary systems is to acquire detailed information about individual planets. With radial 
velocity measurements alone, such information is limited only to the minimum mass Mp sin i 
of the planet, and the semi-major axis a and eccentricity of its orbit. However, if the planet 
also transits its parent star, then it is possible to infer considerably more information. A 
basic transit measurement allows one to infer the radius, mass, and density of the planet, as 
has been done with the only known transiting extrasolar planet, HD209458b (Charbonncau 
et al 2000; Henry et al. 2000). This in turn allows one to place constraints on the planet's 
orbital migration history (Burrows et al. 2000). More detailed photometric and spectroscopic 
data during (and outside of) the transit can be used to study the composition of, and 
physical processes in, the planetary atmosphere (Seager & Sassclov 2000; Seager, Whitney, 
& Sasselov 2000; Charbonncau et al. 2002; Brown, Libbrecht, & Charbonncau 2002), measure 
the oblateness, and thus constrain the rotation rate, of the planet (Hui & Seager 2002; Seager 
& Hui 2002), and to search for rings and satellites associated with the planet (Sartoretti & 
Schneider 1999; Schneider 1999; Brown et al. 2001). 

The 'classical' method of searching for planets via microlcnsing was first proposed by 
Mao & Paczyhski (1991), and subsequently further developed by Gould & Loeb (1992). In 
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this method, a planetary companion to the primary lens star produces a small perturbation 
atop the smooth, symmetric lensing light curve created by the primary. The microlensing 
method has several important advantages over other methods, as well as several disadvan- 
tages (see Gaudi 2003 for a review). The most important advantage is that the strength 
of the planet's signal depends weakly on the planet /primary mass ratio and thus it is the 
only currently feasible method to detect Earth-mass planets (Bennett & Rhie 1996). The 
other advantage is that it enables one to detect planets located at large distances of up to 
several tens of kiloparsecs. However, it also has disadvantages, the most important of which 
is that the only useful information one can obtain is the mass ratio between the planet and 
the primary. Thus classical microlensing searches only allow one to identify the existence of 
the planet, and build statistics about the types of planetary systems, but cannot be used 
to obtain detailed information about the discovered planets. This is especially problematic 
in light of the fact that follow-up of the discovered systems will generally be difficult or 
impossible. 

Recently, Graff & Gaudi (2000) and Lewis & Ibata (2000) proposed a novel method 
of detecting planets via microlensing. They suggested that one could detect close-in giant 
planets orbiting the source stars of caustic-crossing binary-lens events via accurate and de- 
tailed photometry of the binary-lens light curve. In this method, the planet can be detected 
because the light from the planet is sufficiently magnified during the caustic crossing to 
produce a noticeable deviation to the lensing light curve of the primary. The magnitude of 
the deviation is 5^ ~ CpPp where is the ratio of the (unlensed) flux from the planet to 
the (unlensed) flux from the star, and Pp is the angular radius of the planet in units of the 
angular Einstein ring radius of the lens system. The Einstein ring radius is related to the 
physical parameters of the lens system by 



where i?sch = IGMjc? is the Schwarzschild radius of the lens, M is the total mass of the 
lens, D = DogDoi/Dis, and D^s, -Dob As ^-re the distances between the observer-source, 
observer-lens, and lens-source, respectively. For searches in the optical, the hght from the 
planet will be dominated by the reflected light from the star, and ep ~ 10~^ for close-in 
planets.^ Adopting typical parameters, pp ~ 10^'*, and thus Sp ~ 1%. Graff & Gaudi (2000) 
demonstrated that this level of photometric precision is currently within reach of the largest 
aperture telescopes. The exquisite resolution afforded by caustics may allow one to study 



^Because the fraction of reflected light decreases as Cp ex a~^, optical searches will generally only be 
sensitive to close-in planets. However, planets may have significant intrinsic flux in the infrared, enabling 
the detection of more distant companions at longer wavelengths. 
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features on and around the source in detail, and with larger aperture telescopes, one may able 
to study spots and bands on the surfaces of detected planets by looking for small deviations 
to the nominal light curve (Graff & Gaudi 2000). Here we study the signatures of these and 
other structures on lensing light curves, quantify the magnitude of the deviations, and assess 
their detectability using current and future instrumentation. Specifically, we consider the 
signatures of rings and satellites, as well as atmospheric features such as spots, zonal bands, 
and scattering. Lewis & Ibata (2000) considered using variations in the polarization during 
the planetary caustic crossing to probe the composition of the planetary atmosphere. The 
effects of the phase of the planet on the light curve were considered previously by Ashton & 
Lewis (2001). 

The layout of the paper is as follows. In §2, we discuss binary lenses and their associ- 
ated caustic structures, and describe the magnification patterns near caustics. We discuss 
expectations for the existence, stability, and properties of rings, satellites, and atmospheric 
features of close-in extrasolar planets in §3. In §4, we layout the formalism for calculating 
microlensing light curves, and apply this formalism to make a quantitative predictions for the 
deviations caused by planets (§4.1), satellites (§4.2), rings (§4.3), and atmospheric features 
(§4.4). We address the detectability of these deviations in §5, and summarize and conclude 
in §6. 



2. Bineiry Lenses and Caustic Crossings 

If a microlensing event is caused by a lens system composed of two masses, the resulting 
light curve can differ dramatically from the symmetric curve due to a single lens event. The 
main new feature of binary lens systems is the formation of caustics. Caustics are the set of 
positions in the source plane rj) on which the magnification of a point source is formally 
infinite. The set of caustics form closed curves, which are composed of multiple concave 
fine segments that meet at points. The concave segments are referred to as fold caustics, 
whereas the points are cusps. The number and shape of caustic curves varies depending 
on the separation and the mass ratio between the two lens components. Figure 1 shows an 
example caustic structure of a binary lens system with equal mass components separated by 
6e- For more details on the caustic structure of binary lenses, see Schneider & Weiss (1986) 
and Erdl & Schneider (1993). The caustic cross section generally decreases with decreasing 
mass ratio, and decreases for widely and closely separated components. Therefore, the 
majority of caustic-crossing binary-lens events will have caustic structures similar to that 
shown in Figure 1. Most source trajectories (straight lines through the source plane) will 
not intersect the caustic near cusps, therefore the majority of caustic crossings will be simple 
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fold caustic crossings. Near a fold caustic, the total magnification A of a point source is 
generically given by (Schneider, Ehlers & Falco 1992; Gaudi & Fetters 2002), 



where u± is the angular normal distance of the source from the fold in units of 6e, is 
related to the local derivatives of the lens potential, and describes the effective 'strength' of 
the caustic (also in units of ^e), and Aq is the magnification of all the images unrelated to 
the caustic. The divergent nature of the magnification (as A oc u]^^"^) translates into high 
angular resolution in the source plane. 

Due to the divergent magnification near a caustic, the light curve of an caustic-crossing 
binary lens event is characterized by sharp spikes which are generally easily detectable. Mao 
& Paczyiiski (1991) predicted that ~ 7% of all events seen toward the Galactic bulge should 
be caustic-crossing events. Analyses of the databases of the MACHO (Alcock et al. 2000) and 
OGLE ( Jaroszynski 2002) lensing surveys have demonstrated that binary microlensing events 
are being detected at a rate roughly consistent with theoretical predictions. Therefore a sig- 
nificant sample of caustic-crossing binary-lens events should be available each year for planet 
searches. The planetary caustic crossing is expected to occur within ~ 1 day (a/0.1 AU) of 
the stellar caustic crossing. Therefore, in order to detect close-in planets and their associated 
structures, extremely dense sampling for a period of ~ 2 days centered around one of the 
stellar caustic crossings is required. Caustic crossings occur in pairs, and although the first 
caustic crossing may not be detected real time because of its short duration, it can be inferred 
afterwards from the enhanced magnification interior to the caustic. Followup observations 
can therefore be prepared before the second caustic crossing, and dense sampling throughout 
the second caustic crossing will be possible. 

The usefulness of caustic-crossing binary-lens events has already been demonstrated in 
numerous ways. Precise photometry during the caustic crossings of several events has been 
used to measure the limb darkening profiles of stars located in both the Galactic bulge and 
the Small Magellanic Cloud (Albrow et al. 1999; Afonso et al. 2000; Albrow et al. 2001; An 
et al. 2002), and spectra taken during the unusually long caustic crossing of one event has 
been used to resolve the atmosphere of a X-giant in the bulge (Castro et al. 2001; Albrow 
et al. 2001). It also has been proposed that irregular structures on the source star surface 
such as spots can be studied in detail by analyzing the light curves of caustic-crossing binary 
lens events (Han et al. 2000; Chang & Han 2002). The same principles that make these 
measurements possible also allow one to study tiny structures on and around the planetary 
companions to the source stars of caustic-crossing events. 
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What is the ultimate resolution that can be obtained during a caustic crossing? In 



-6- 



geometric optics, this is set by the frequency of measurements and the number of photons 
that can be acquired in a given measurement, so that extremely small structures in the 
source plane could, in principle, be probed with sufficiently large telescopes and sufficiently 
high cadence. However, for very small sources, geometric optics breaks down, and diffraction 
effects become important. For a fold caustic crossing, this occurs when the angular size 
of the source in units of 6'e satisfies ps < Pdi where (Ulmer & Goodman 1995; Jaroszynski & 
Paczynski 1995), 



Here A is the wavelength of the hght. For A ~ 800 nm (7-band) and M — O.SM©, ~ 
7 X lQ~^Ur^^^. For nearly equal-mass binaries, such as shown in Figure 1, Ur ~ 1. Thus 
the ultimate resolution achievable during a caustic-crossing is p^O^ ~ 20 pas, or a length of 
PdO^Dos ~ 26 km at the distance of the bulge. This is considerably smaller than the size of 
any of the structures we consider here, so we can safely ignore diffraction effects. 

3. Satellites, Rings, and Atmospheric Features of Close-in Planets 

All of the planets in our solar system have at least one satellite, with the exceptions of 
Mercury and Venus, so it seems at least plausible that satellites are common by-products of 
the formation of planetary systems. Satellites are perturbed by the tidal bulges they induce 
on their parent planets, and their orbits evolve under the influence of this torque. If the 
timcscale for this evolution is sufficiently short, the satellite will either spiral inward until 
impacting with the parent planet, or outward until it reaches the Hill radius of the planet, and 
is lost to the parent star. The survival of satellites has been considered by numerous authors 
in the context of our solar system (see, e.g. Ward & Reid 1973). More recently, Barnes 
& O'Brien (2002) studied the lifetimes of sateUites in the context of extrasolar planetary 
systems, and found that satellites with Mgat > -^e cannot survive for more than ~ 5 Gyr 
around Jupiter-like planets with separations a < 0.25 AU, assuming a solar-mass primary 
star. Therefore relatively massive satellites around close-in planets are expected to be rare. 
However, the lifetime of satellites depends sensitively on the (uncertain) tidal Q value of the 
planet, and therefore may be significantly in error. Furthermore, planets around lower-mass 
stars are expected to retain their satellites for a longer period of time. 

Satellites of close-in planets are unlikely to have substantial atmospheres because their 
surface escape speeds arc generally small enough that most light element gases will have 
evaporated over the age of the system. The equilibrium blackbody temperatures of satellites 
of planets with a < 0.1 AU are T > 900 K, and any gases with an atomic mass less than 
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~ 18 (including H, He, CH4, NH3 and H2O), will have been entirely lost over the ~ 5 Gyr 
age of the system for satellites with mass < M®. Therefore, the scattering surface of any 
extant satellite will likely be rocky. 

Ring systems exist around all of the gas giants in our solar system, and therefore might 
also be expected to be common debris from planetary formation. The rings we observe 
in the solar system have varied properties, but at least some, such as the rings of Saturn, 
are composed of icy materials which give rise to high albedos. Such high albedos would 
aid considerably in detection in the current context. Unfortunately ices cannot exist as 
separations less than. 



where is luminosity of the parent star and Tsub — 170K is the sublimation temperature 
of ice. Thus close-in planets cannot have rings composed of icy material. Rocky rings are 
not precluded; however the constituent particles are subject to numerous dynamical forces, 
including Poynting-Robertson (PR) drag, viscous drag from the planet exosphere, torques 
from satellites and/or shepherd moons, and internal collisions. The characteristic decay time 
for PR-drag is (Goldreich & Tremaine 1982) 



where p and r are the density and radius of the particles, respectively. For close-in planets, 
the decay time ty^ for viscous drag is likely to be considerably larger than (Goldreich & 
Tremaine 1982), unless the planet's exospheres are quite dense, p > 10~^^ g cm~^. Interpar- 
ticle forces only serve to spread the ring. Therefore, the dominant effect, aside from satellite 
perturbations, is likely to be PR drag. It is clear that rings of close-in planets will be lost 
to the planet on a relatively short timescale unless they are stabilized by interactions with 
satellites. Since, as we have just discussed, satellites around close-in planets are themselves 
generally not long-lived, it is not clear that this is a viable method of maintaining rings. A 
definitive exploration of the stabihty of rings and sateUite systems around close-in planets is 
beyond the scope of this paper, but warrants future study. 

An enormous amount of effort has been devoted to modeling of the atmospheres of 
extrasolar planets, with ever increasing levels of sophistication (see, e.g. Saumon et al. 1996; 
Burrows et al. 1997; Seager & Sasselov 2000). Special emphasis has been placed on close- 
in planets (Seager & Sasselov 1998; Goukenleuque et al. 2000). The recent confrontation of 
observations of the radius and atmosphere of HD209458b (Brown et al. 2001; Charbonneau et 
al 2000) with theoretical predictions (Seager & Sasselov 2000; Burrows et al. 2000) generally 
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indicates that considerably more work needs to be done (Guillot & Showman 2002; Fortney 
ct al. 2003). The problem of modeling the atmospheres of close-in planets is an especially 
difficult one, due to the fact that these planets are tidally locked to their parent stars, and 
subject to strong stellar irradiation. It seems likely that non-equilibrium processes, weather, 
and photoionization will all play a role in accurate theoretical models. It is therefore perhaps 
a bit premature to speculate on the existence and nature of surface features, such as zonal 
bands and spots, on extrasolar planets. There are intriguing indications, however, that close- 
in planets may possess large-scale surface features. Showman & Guillot (2002) and Cho et al. 
(2003) showed that the extreme day-night temperature difference in close-in, synchronized 
planets drive large-scale zonal winds that can reach ~ 2 km s~^. These circulation patterns 
result in non-uniform surface temperatures which may lead to significant variations in the 
scattering, refiecting, and absorbing properties of the atmosphere. 

Regardless of whether or not surface irregularities exist on extrasolar planets, it is 
clear that a uniform surface brightness, which has been assumed in previous microlensing 
simulations (Graff & Gaudi 2000; Lewis & Ibata 2000; Ashton & Lewis 2001), will likely 
not be an accurate representation of the global illumination pattern of the planet. Even for 
simple Lambert scattering, in which each area element of the planetary atmosphere reflects 
the incident flux uniformly back into the 27r available solid angle, the surface brightness profile 
of the planet is non-uniform due to projection effects. Departures from Lambert scattering 
are expected, and depend on such properties as the particle size of the condensates in the 
atmosphere (Seager & Sasselov 2000). 

To summarize, it remains unclear whether satellites or ring systems can exist around 
close-in extrasolar planets. Models of the close-in planetary atmospheres have not yet reached 
the level of sophistication required to definitively predict whether large-scale surface features 
such as zonal bands or spots will be present. It seems quite likely, however, that the surface 
brightness profile of extrasolar planets will not be uniform, as has previously been assumed 
when calculating the effects of microlensing. We will therefore proceed with rampant opti- 
mism, and assume that all of the above structures, (rings, moons, zonal bands, spots, and 
non-uniform surface brightness profiles) may exist, and consider the nature and magnitude 
of their effects on microlensing light curves. 



4. Quantitative Estimates 

In this section, we estimate the magnitude of the features in lensing light curves produced 
by satellites, rings, and atmospheric features, relative to the nominal light curve produced 
by an isolated, circular planet with uniform surface brightness. For the most part, we use 
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semi-analytic means to produce quantitative estimates in order to elucidate the dependence 
of the deviations on the parameters of the system. However, in some cases the signals cannot 
be computed analytically. We therefore complement our semi-analytic results with detailed 
numerical simulations. 

The time-dependent total flux F{t) from a system composed of a star and N additional 
source components being microlensed can be generically written as, 

N 

F{t) = F,A,{t) + J2 FiMt) + B, (6) 

i 

where is the unlensed flux of the star, A^, is the magnification of the star, Fj and Ai are the 
unlensed flux and magnification of the iih. additional component (which may include planets, 
satellites, rings, etc.), and B is any unlensed blended fiux. We will henceforth assume no 
blending, but note that any deviation in the light curve is suppressed in the presence of 
blending by a factor (1 + fb/Ani,)^^, where = B/{F^ + ^ F,) is the ratio of blend to total 
source flux, and Anb is the total magnification in the absence of blending. 

Generally, wc have that Fi ^ F^,, and <C F*, and thus the observed magnification 

can be approximated as 

F(t) 

Aobs = p , V; p ^A, + AAtot; AAtot = J2 (7) 

where we have defined AAtot, the extra magnification due to the planetary companion and 
associated structures, and the fiux ratio fi = Fi/F^ between the ith component and the star. 

For a uniform, circular source sufficiently close to a hnear fold caustic, the magnification 
is (Chang 1984; Schneider, Ehlers & Falco 1992), 

Af,^(^^yGo(z)+Ao, (8) 

where z = uj_/ ps is the normal distance of the source from the caustic in units of the source 
size. The function Go{z) describes the normalized light curve for a uniform, circular source 
crossing a fold caustic, and can be expressed in terms of elliptic integrals (Schneider & Weiss 
1986). It is useful because it can be used to describe the magnification of any source that 
can decomposed into components with azimuthal symmetry. Gq has a maximum of ~ 1.38 
at z ~ —0.66. In the range — 1 < ^; < 1, the mean value of Go{z) is 0.95 and the RMS is 
1.04. The magnification can also be written as a function of time by defining tr = UrtE esc 7, 
and z — {t — tec)/ At, where tec is the time when the center of the source crosses the caustic, 
and 7 is the angle of the trajectory with respect to the caustic, and the timescale of the 
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caustic crossing is At = ps^ECSC7. Here — 9^/ ^ is the timescale of the primary event, 
and is the relative lens-source proper motion. 

For typical microlensing events toward the Galactic bulge, — 320//as and ji ~ 
12.5 km s~^kpc~^, and thus t-p. ^ 22 days. We will assume that the primary source is a 
G-dwarf in the bulge, i.e. that it has a radius — Rq and is located at the distance 
Dos = 8 kpc. The angular radius is 9^ ~ 0.6//as, and thus the dimensionless source size is 
p* — 9^/9e — 1.8 X 10~^. The caustic crossing timescale for a source of physical radius R is, 

AU ^ 1 hr CSC 7 ^-^^ . (9) 
Thus the primary caustic crossing is expected to last 2At^ ~ 2 hr for R^ — Rq. 



4.1. Planet 

The largest contribution to AAtot generally will be from the planet itself, AAp = fpAp. 
In the case of light reflected by a planetary atmosphere or surface, the flux fraction /p will 
depend on the radius Rp of the planet, its distance a from the star, the scattering properties 
of the atmosphere, and the phase of the planet. Generically, the flux ratio /p between the 
planet and star can be written as (Sobolev 1975), 

/p = ep$(a), (10) 

where a is the phase angle, defined as the angle between the star and Earth as seen from 
the planet, $(q;) is the phase function, and Cp is the fiux ratio at a = 0, 

e,^^,(^)^2.28xlO-A(^)-^(|H_)= (11) 

Here Ap is the geometric albedo of the planet. For a Lambert sphere, Ap — 2/3, and 

1 

$(«) = — [sin a + {% — a) cos a] . (12) 

TT 

The magnification Ap of the planet will depend on the size of the planet, as well as on 
its surface brightness. The surface brightness, in turn, depends on the phase of the planet, 
as well as the scattering properties of the atmosphere. The effects of the phase of the planet 
on the magnification have been considered by Ashton & Lewis (2001), and we consider the 
effects of Lambert-sphere scattering on the magnification in §4.4.3. We will therefore assume 
that the planet has a uniform surface brightness and is at full phase {a — and thus ^ — 1), 



unless otherwise stated. In this case, the magnification of the planet is simply given by 
equation (8), with ps = pp, where Pp is the angular radius of the planet in units of ^e- 
Adopting this expression, the contribution of the perturbation from the planet is, 

AAp^epi^-^j Goizp). (13) 

Note that we have assumed that fpAo <^ 1. In most cases, Aq is of order unity, and this will 
be an excellent approximation. We will furthermore assume that JiAq -C 1 in deriving all 
analytic expressions. During the planetary caustic crossing, the RMS^ of Gq is ~ 1. Thus 
the magnitude of the deviation is given by the coefficients of the Go function in equation 
(13), ^ 

For binary lenses with caustic structures similar to that shown in Figure 1, Ur is of order 
unity (Lee, Chang, & Kim 1998). 

Thus for typical microlensing bulge parameters, Sp — 1.67x10'^ Ap{Rp/ Rjupj^^^ia/O.l AU) 
and for a planet with properties similar to HD209458b {Rp — l.SATRjup and a — 0.0468 AU, 
Brown et al. 2001), and = 2/3, Sp ~ 0.8%. From equation (9), the duration of planetary 
caustic crossing is 2Atp ~ 16 min esc 7 for Rp — 1.34:7 Rjup- 



4.2. Satellites 

The magnitude of the deviation caused by a satellite can be estimated using the same 
formalism as used for the planet (§4.1). As for the planet, we will assume that the sateUite 
has a uniform surface brightness and phase a — 0. The deviation AAgg, caused by the satellite 
is then, 

/ \ 1/2 

A^a^esa — Go(0, (15) 

VPsa/ 

where is the fiux ratio between the satellite and star, and is the dimensionless size of 
the satellite. In analogy to the case of the planet alone, we can define 5sa = £sa{ur/ PsaY^'^ to 
be the magnitude of the deviation from the satellite. If we assume that the distance of the 
satellite from the planet is small compared to a, we can relate this to the magnitude of the 



^The RMS is the relevant quantity for signal-to-noise considerations, see §5. 
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deviation due to the planet by, 




) 



3/2 



'sa 



(16) 



where Kgs, = Asa./Ap is the ratio of the geometric albedos of the satellite and planet, and i?sa 
is the radius of the satellite. The ratio will depend quite strongly on the compositions of 
the planet and satellite. As we argued in §3, satellites of close-in planets are unlikely to have 
substantial atmospheres, and therefore the scattering surface will be rocky in composition. 
Albedos of rocky bodies depend on their composition, but are generally low, < 0.3. For 
definiteness, we will assume that Kga = 0.1, and assume a satellite radius of Rsa, — 0.2Rp. 
This is ~ 2i?0 for a Jupiter-size planet. Then Sga. — l%Sp. 

Assuming typical bulge parameters, an analog of the HD209458 system at D^g = Skpc, 
and for the second caustic crossing of the dashed trajectory in Figure 1, which has the 
properties 7 = 71° and Ur = 0.82, we find Sp = 7.19 x 10~^ and Atp = 8.48 min. Adopting 
these parameters. Figure 2 shows the total magnification associated with the planet and 
satellite, AA^ot — AAp + AA^a,, as well as the extra magnification from the satellite alone, 
AAsa- We show the effect for various satellite positions. The time it takes for the caustic to 
cross the satellite is 2Atsa — 3.4 min. 



Rings of extrasolar planets have such low mass that they have no observable dynamical 
effect on the host star's motion. However, as shown by the example of Saturn's ring, they can 
be significantly more extended than planets. This makes them much easier to be identified 
by transits, for which the signal varies as the area Q of the feature that occults the star, and 
microlensing, for which the signal is oc due to the competing effects of the amount of 
refiected light (oc fi), and the magnification (oc f2~^/''). 

We can obtain an analytic estimate of the signal caused by a ring by assuming a face-on 
geometry. We model the ring as a circular annulus with uniform surface brightness, outer 
radius Rout, and inner radius The magnification from the ring is then. 



4.3. Rings 



AA,i = 5,iH, 



(17) 



where we have defined 




3/2 



3/2 



H = 



(18) 
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and Zout = Zp{Rp/ Rout), and similarly for Zm- The shape and magnitude of the function H 
will depend on the relative sizes of Rout, Rin, and -Rp, but will generally be of order unity for 
the ring systems considered here. The magnitude of the deviation due to the ring is therefore 
roughly S^, which in terms of the deviation from the planet is, 

Sri = KriSp, (19) 

where K^-i = An/ Ap is the ratio of the geometric albedos of the ring and planet. We argued 
in §3 that any surviving ring systems of close-in planets must be rocky in nature, and thus 
the albedos will generally be small. Adopting n^i = 0.15, S^i = 15%6p. 

Figure 3 shows the total magnification associated with the planet and ring system, 
AAfot = AAp+AAj.[, where wc have adopted the same parameters for the planetary deviation 
as in §4.2, namely 6p = 7.19 x 10~'^ and Atp = 8.48 min. We also show the magnification 
from the ring system alone, AAri. We show the effect of the ring for various inner radii, and 
in order to isolate the effect of varying the inner radius on the resulting light curve, we fix 
the total area of the ring systems to be ^{Rl^^ — R^^) — 10i?p. The two small bumps on 
the left and right sides of the primary peak are caused by the ring's entrance and exit of 
the caustic. The time between the two ring-induced bumps (or between one of the bumps 
and the primary peak), relative to the time scale of the planetary perturbation, is a measure 
of the relative dimension of the ring compared to the size of the planet disk. The time it 
takes for the caustic to cross the ring system is 2Atr\ = 2{Rout/ Rp)Atp. For the largest ring 
system shown in Figure 3, this is ~ 74 min. We find that the signal of the ring generally 
decreases slowly as the gap between the planet and the ring increases. 

Although our analytical expressions for the deviations induced by a simple, face-on ring 
system are useful in that they allow one to gain insight and find relatively simply scalings 
for the magnitude of the effect, they are limited in their scope. In particular, they cannot be 
used to access the effects of inclined ring systems. For this, numerical integration must be 
employed in order to calculate the magnification of the ring. We consider sources crossing 
the second caustic of the dashed trajectory in Figure 1. To calculate the light curve, we first 
compute the full binary-lens magnification on each area element on the surface of the source 
and then average the magnifications of the individual elements, weighting by the surface 
brightness of each element. Wc again assume typical microlensing parameters, and a source 
system analogous to HD 209458 at 8 kpc. This yields Cp = 1.26 x 10^"^, and ppt^ = 8.02 min. 
The effects of Ur and 7 are included implicitly in our numerical integration of the binary- 
lens magnification for the specific trajectory we have adopted. For the planet-only case, we 
find our numerical and analytic light curves agree quite well, indicating that our numerical 
integrations are accurate, the caustic is well- approximated by a simple linear fold, and that 
we are using the appropriate values of 7 and Ur in our analytic expressions. 
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In the numerical simulations, we model ringed planets in the same manner as the analytic 
calculations. We assume that the ring is an infinitely thin annulus without any gap, with 
inner and outer radii i?j„ and Rout, respectively. Wc assume the planet has a uniform albedo 
of Ap = 2/3 for pure Lambert scattering. We assume that the ring has an albedo of A^i = 0.1, 
and thus = 0.15, unless otherwise specified. For the simulations, we take the effects of 
the planet's phase and the ring's inclination into consideration. For some specific geometries 
of the planet, host star, and the observer, the planet will cast a shadow on the ring. We 
also take this effect into consideration. We then investigate the variations of the pattern of 
ring-induced anomalies depending on these various factors affecting the shape of the ring. 

Defining the projected shape of a ringed planet requires many parameters, such as the 
inclination of the ring, the phase angle, the radius of the planet disk, and the inner and 
outer ring radii. As a result of the large number of parameters, it is often difficult to imagine 
the planet's shape based on these parameter values. We therefore simply use small icons to 
characterize the planet and ring shapes instead of specifying all these parameters whenever 
we present light curves resulting from specific realizations. 

Figure 4 shows the effect of the width of the ring on the light curve. All of the systems 
have a common inner ring radius and gap between the planet disk and the ring, but different 
outer ring radii. We have also assumed that the ring system is viewed at an angle of i = 75° 
with respect to the normal to the ring plane. Not surprisingly, increasing the width of the 
ring increases the magnitude and duration of the ring-induced perturbation. 

Figure 5 shows the dependence of the ring-induced anomaly pattern on the albedo of 
the ring particles. We test three different albedos of = 0.05, 0.1, and 0.2 (corresponding 
to = 0.075, 0.3, and 0.6), and the difference in the greyscale of the rings in the icons 
represents the variation of the albedo. As expected, the magnitude of the ring signal is 
proportional to the albedo of ring particles. 

In Figure 6, we demonstrate the effects of a shadow cast on the ring by the planet. In 
this case, we must also take into account the phase of the planet to be self-consistent. For 
the geometry depicted in Figure 6, the planet is at quarter phase. We find that the shape 
of the part of the light curve arising from the planet is strongly dependent on the phase of 
the planet, as discussed in detail by Ashton & Lewis (2001), however the shape of the signal 
due the ring does not depend strongly on the effect of the shadow of the planet, due to the 
relatively small surface area of the ring occulted by the planet. 

From the light curves in Figures 4-6, we conclude that the typical magnitude of 
ring-induced deviations is S^i ~ Uri^p, confirming our analytic estimates, although there are 
considerable variations depending on the inclination and size of the ring. 
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4.4. Atmospheric Features 

In contrast to the signatures of satellites and rings, atmospheric features that are located 
on the surface of the planet will only produce detectable deviations while the planet is 
resolved. The planet is only effectively resolved when it is within about one planet radius 
from the caustic. Therefore, deviations caused by spots, bands, or otherwise non-uniform 
surface brightness profiles will only be noticeable during a time ~ 2Aip centered on the 
caustic crossing; outside of this the flux from the planet will essentially be given by the 
unresolved fliix (i.e. the mean surface brightness times the area of the planet disk) times the 
magnification of a point-source at the center of the planet disk. 



4.4- 1- Spots 

Spots, such as the Great Red Spot on Jupiter, have been observed on the gas giants in 
our solar system, and are regions of cyclonic activity that have slightly different temperatures 
and pressures from their surrounding atmospheres. As a result, the colors and albedos of 
spots are also slightly different. These spots can be quite large relative to the planetary radii; 
the Great Red Spot has a size larger than R^. In this subsection, we provide an analytic 
estimate of the effect of a spot on the microlensing light curve. We model the spot as a circle 
of radius i?sp, with an albedo equal to a fraction of the albedo of the remainder of the 
planetary surface. The deviation caused by the spot is then. 



AA. 



sp 



sp 



1 - (1 - «sp) {R.jRpf 



^)'Go(.p)-(^y^'Go(.sp) 



(20) 



Note that, in deriving equation (20), the mean surface brightness of the spotted planet has 
been normalized to that of the planet without the spot. This ensures that the magnifications 
of the two cases are identical when the source is not resolved. That this is true can be seen by 
noting that, for ^^gp, Zp — > — oo, the term in square brackets goes to zero, because Go{zsp) — > 
{Rsp/Rpy^'^Go{zp). Since we are generally concerned with small spots with Rsp/Rp <^ 1, we 
can make an estimate of the magnitude Ssp of the deviation caused by the spot by ignoring 
terms in equation (20) of order {Rgp/RpY or higher, and looking at the resulting coefficient 
to the Gq(zsp) function. We find that 



Ssp ~ -(1 - r,p) ( ^ ) dp. 



Rr 



3/2 



(21) 



For Ksp = 0.8, and Rsp = 0.2Rp (~ 2R^ for Rp = Rjup), we find 5sp ^ 2%5p. 
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In Figure 7, we show the total magnification AAtot, with and without a spot of radius 
Rsp = 0.2Rp and relative albedo Kgp = 0.8. We have adopted the same parameters as in 
§§ 4.2 and 4.3, 6p = 7.19 x 10"^ and Atp = 8.48 min. We vary the position of the spot 
as shown. We also show A Asp, the deviation caused by the spot from the light curve of a 
planet with uniform surface brightness. We find the magnitude of the deviation caused by 
the spot to be ~ l%5p, in rough agreement with our analytic estimate. Since, for linear 
fold caustics, the magnification is independent of position along the direction parallel to the 
caustic, the results in Figure 7 are applicable for any circular spot located on the star at the 
same perpendicular distance from the caustic as the spots shown. Microlensing of spotted 
planets is analogous to microlensing of spotted stars, see Han et al. (2000), Lewis (2001), 
and Chang & Han (2002) for examples of such hghtcurves. 



4- 4- 2- Zonal Bands 

As can be seen on the surface of Jupiter, gaseous giant planets may exhibit color varia- 
tions on their surface, e.g. zonal bands, which will cause surface brightness variations within 
a given spectral band. In this subsection, we investigate whether zonal bands can produce 
noticeable signatures in lensing light curves. 

We model zonal bands by stripes which are parallel with the equator of the planet. We 
assume that the albedos of these stripes alternate with relative values K^g,, and we vary the 
total number of stripes. Analytic estimates are generally impossible for arbitrary inclinations 
of the planet, however, for the special case when the planet is seen pole-in, the pattern of 
the zonal bands is simply concentric annuli with alternating albedos. In this case, we can 
find a semi-analytic expression for the deviation from a uniform surface brightness due to 
the zonal bands. The resulting expression is somewhat complicated. 



1 - «ba 



l-(l-«ba)(i?ba/i?p)^ 



-H' 



(22) 



where 



H' = 



R 



ba 



Rn 



A^ba-l / p N 3/2 

GoiZp) - E (-1)^+^ [^j Go(.ba.), 



(23) 



and A^ba is the number of bands, i?ba,t is the outer radius of the ith band, z^^i = ^p(-Rp/-Rba,i)) 
£ = 1 if A'ba is odd, and £ = if Aba is even, and 



i=l 



2 

ba.i" 



(24) 
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Note that, as A^ba oo, A^ba 0. The term in square brackets in equation (22) cannot be 
reduced to a simpler expression, and must be calculated explicitly using the analytic form 
for the Go function. 

In Figure 8, we show the total magnification AAtot of the planet, with and without A'ba 
bands with relative albedos Kba = 0.8, similar to that of the Jupiter (Pilcher & McCord 1971). 
We have adopted the same parameters as in the previous sections, except here we consider 
the second caustic crossing of the solid trajectory in Figure 1, which has the properties 
Ur = 0.31 and 7 = 89°. This yields 6p = 4.41 x 10"^ and Atp = 8.02 min. We vary the 
number of bands from iVba = 3 — 5, and the bands have equally-spaced radii (and thus cover 
different surface areas on the source). We also show AAba, the normalized deviation caused 
by the banded structure from the light curve of a planet with uniform surface brightness. 
We find the magnitude of the deviation caused by the bands to be 

5ba ~ 30%(1 - «:ba)A^b?<^P- (25) 

Here the scaling with ^ba is only approximate. The scaling with A^ba depends on the geometry 
of the zones; P — 1 ior zones with equally-spaced i?ba,i (as shown in Fig. 8), whereas (3 — 1/2 
for equal-area bands. For equal-area zones, the numerical coefficient in equation (25) is also 
somewhat smaller, ~ 20%. 

For geometries where the planet is not pole-on, we must resort to numerical calculations. 
For these calculations, we assume a total of nine bands (with four dark lanes), with relative 
albedos of Kba = 0.8, as in the previous example. The albedos are normalized such that the 
average albedo is 2/3. Figure 9 shows the effect of zonal bands for a planet with inclination 
i = 90 (i.e. the axis of rotation in the plane of the sky), and various orientations of the 
equator with respect to the caustic. As before, we assume the source trajectory indicated by 
the solid line in Figure 1. The solid curve is the light curve resulting from a planet having a 
uniform surface brightness with ^p. From the figure, one finds that the deviations induced 
by the zonal bands are ~ l%5p, similar to the pole-on case. Note that these deviations are 
generally an order of magnitude smaller than the typical deviations induced by rings. 



^.^.5. Lambert Sphere Scattering 

Even without any irregular structure, the surface brightness profile of a planet will 
generally not be uniform due to projection effects, and the scattering properties of the 
atmosphere. The surface brightness distribution will therefore vary depending on the latitude 
■0 and longitude u of the planet's surface as well as the planet's phase angle, a. To illustrate 
the pattern of lensing light curve deviations caused by a realistic atmosphere, we adopt the 
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simple assumption of pure Lambert scattering, where the incident radiation from the host 
star is scattered isotropically. Under this assumption, it can be shown that the surface 
brightness profile is 

S{a,u,^) = S- , S=—, (26) 

for u > a — 7i/2 and zero otherwise (Sobolev 1975). Here S is the mean surface brightness 
of the planet, and $(a) is the phase integral introduced in §4.1 and displayed explicitly in 
equation (12). We show the surface brightness distribution for a Lambert sphere with a — 
in Figure 10, along with the resulting microlensing light curve. As in §4.4.2, we have assumed 
the solid source trajectory in Figure 1. We compare this to the light curve resulting from a 
planet with a uniform surface brightness and an albedo equal to the geometric albedo of the 
Lambert sphere (i.e. Ap = 2/3). We find that the light curve of a Lambert sphere differs from 
a uniform surface brightness profile by ~ 10%5p. The light curve from a Lambert sphere is 
more highly magnified, due to the fact that the surface brightness profile is more centrally 
concentrated, and therefore the source is effectively smaller. The precise shape of the light 
curve from the planet will depend on the scattering properties of the atmosphere, which 
in turn depend on the constituents of the atmosphere, such as the size of the condensates 
(Seager & Sasselov 2000). Therefore resolution of the planetary caustic crossing would 
provide invaluable information about the physical processes in the planetary atmosphere. 
This is complementary to the suggestion by Lewis & Ibata (2000) of probing the planetary 
atmosphere via polizarization monitoring during the planetary caustic crossing. 



5. Detectability 

In this section, we review the magnitudes of the effects of the features we have considered, 
and assess their detectability with current and/or future telescopes. Table 1 summarizes our 
analytic expressions for the magnitudes 5x of the deviations caused by each structure x that 
we considered (satellites, rings, spots, and zonal bands), in terms of the magnitude Sp of 
the deviation due to only the planet. Also shown is the characteristic timescale At^ of each 
deviation, relative to the timescale of the planetary caustic crossing Atp. 

By approximating the perturbations from each structure x as boxcars with amplitudes 
Sx and durations At^, we can write down approximate expressions for the ratio of the signal- 
to-noise Qx for a given perturbation to the signal-to-noise of the planetary pertrubation. 



X 



At 



Atp 



1/2 



(27) 
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These ratios are displayed in Table 1; they allow one to easily estimate the detectability of 
the various features in terms of the detectability of the planetary deviation. For example, 
if one were to require a signal-to-noise ratio of Qp = 20 for a secure detection of the planet 
signal, then signal-to-noise of the deviation from a ring with relative albedo Kj-i = 0.15 and 
relative (outer) radius of Rn/ Rp = 4 would be — 0.15 x ^/A x 20 = 6. For reasonable 
parameters, we expect that Qx/Qp -C 1 for satellites, spots, and bands, whereas for rings 
Qri/Qp - 30%. 

We now provide a more accurate estimate of the expected signal-to-noise for the various 
features. Assume that a microlcnsing light curve is monitored continuously from t,nin to 
tmax, for a total duration T, with a telescope that collects photons per second per unit 
flux. The signal-to- noise ratio Qx of a deviation AAx{t) is then, 

where is the unlensed flux of the primary star, and B is the background flux (sky -|- 
unresolved stars) . The term in curly brackets is essentially the RMS of the deviation during 
the time of the observations. We assume that the source system is an analog of HD20458 
at Dos = 8 kpc. The primary has = 19.5 (a GOV star at 8 kpc with 1.2 magnitudes of 
extinction), and its planet has Rp = 1.347Rjup, a = 0.0468 AU, and Ap = 2/3. Adopting 
typical bulge parameters (M = O.SAf^.Doi = 6 kpc, Dqs = 8 kpc), and a caustic crossing 
with properties Ur = I and 7 = 90°, this gives 6p = 7.93 x lO"'^ and Atp = 8 min. We assume 
a total background flux of 19.3, which includes the moon-averaged sky background at an 
average site, and the contribution expected from unresolved stars in the bulge for a seeing of 
0.75". We assume that a telescope of diameter At collects rii, — 2700(^4^/10 m)^ photons per 
second at / = 20, which corresponds to an overall throughput (including detector efficiency) 
of ~ 50%. Finally, we assume that the light curve is monitored from — 5Atp before the 
caustic crossing until 5 Atp after the caustic crossing of the primary star, for a total duration 
of T ~ 80 min. 

The resulting signal-to- noise values for the various deviations are tabulated in Table 1. 
For lOm-class telescopes, the deviation from the planet should be detectable with Qp ~ 15. 
This is in rough agreement with the results of Graff & Gaudi (2000) and Ashton & Lewis 
(2001) for full phases. However, Ashton & Lewis (2001) find that the signal-to-noisc depends 
quite strongly on the phase of the planet. Adopting a different phase would therefore affect 
the relative signal-to-noise Qx/Qp between the planet and the satellite, ring, spot or band 
features, but generally not the absolute signal-to- noise Qx- For the deviations from the other 
structures, we have adopted the parameter values appropriate to the short-dashed line in 
Figure 2 for the sateUite, the dotted line in Figure 3 for the ring, the short-dashed fine in 
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Feature Magnitude Timescale Relative S/N Absolute S/N 

^t^/^tp Qx/Qp lOm 30m 50m 100m 



Planet 


= (^) 


1/2 


15.1 45.2 75.4 150.7 


Satellite 

Ring 

Spot 
Bands 


/ \ 3/2 / \ 2 

f.. 1 Rsa \ -Rsa 1^ 1 Rsa \ 

^ . Rri *i- . ^ri ^ ^ 
Rp \RpJ 

\3/2 N 2 

0.3(1 - K^,)N-f 1 0.3(1 - K^,)N-f 


0.1 0.2 0.3 0.7 

6.1 18.4 30.7 61.4 

0.1 0.3 0.5 0.9 
0.1 0.3 0.6 1.1 



Table 1 Estimated Signal-to- Noise Ratios for Planetary Structures. 



Figure 7 for the spot, and the long-dashed line in Figure 8 for the zonal bands. For these per- 
turbations we find signal-to-noise ratios of Qsa = 0.1 (satellite), Qri = 6.1 (ring), Qsp = 0.1 
(spot), and Qba = 0.1 (zonal bands). These are in rough agreement with the expected scal- 
ing with Qp, and generally indicate that it will be impossible to detect spots, bands, and 
satellites with lOm-class telescopes. Rings are potentially detectable, but only under some- 
what optimistic scenarios, i.e. large, face-on rings. We therefore consider the detectabihty 
with larger-apertmc telescopes, such as the proposed 30- meter aperture California Extremely 
Large Telescope (CELT) (Nelson 2000), or the European Space Agency's proposed lOO-meter 
aperture Overwhelmingly Large Telescope (OWL) (Dierickx & Gilmozzi 2000). We find that 
rings will generally be detectable with reasonable signal-to-noise for > 30 m, and spots, 
satellites and bands are likely to be undetectable with any foreseeable telescope. 

For the deviations caused by Lambert scattering shown in Figure 10, we find that 
Q ~ l%Qp. Therefore, the non-uniform nature of the surface brightness may be measurable 
with lOOm-class telescopes. 



6. Summary and Conclusion 

Planetary companions to the source stars of caustic crossing microlensing events can 
be detected via the brief deviation created when the caustic transits the planet, magnifying 

—1/2 

the refiected light from the star. The magnitude of the planetary deviation is 6p ~ eppp , 
where Cp is the fraction of the fiux of the star that is refiected by the planet, and Pp is the 
angular size of the planet in units of the angular Einstein ring radius of the lens. For giant, 
close-in planets (similar to HD20958b), ep ~ 10~^, and for typical events toward the Galactic 
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bulge, pp ~ 10 ^. Thus Sp ~ 1%, which is accessible to lOm-class ground-based telescopes. 

Due to the extraordinarily high angular resolution afforded by caustic crossings, fine 
structures in and around the planet are, in principle, also detectable. We first presented a 
brief discussion on the existence and stability of satellites, rings and atmospheric features 
of close-in planets, concluding that although rings and satellites may be short-lived due to 

dynamical forces, the ultimate fate of such structm^cs is not clear. There are good reasons 
to believe that atmospheric features may be important in close-in planets. We therefore 
considered the signatures of satellites, rings, spots, zonal bands, and non-uniform surface 
brightness profiles on the light curves of planetary caustic-crossings. Where possible, we 
used semi-analytic approximations to derive useful expressions for the magnitude of the 
deviations expected for these features, as a function of the relevant parameters, such as the 
albedo or size of the feature. We express these deviations in terms of 5p, the magnitude of 
the planetary deviation. 

We find that rings produce deviations of amplitude ~ 10%Sp, whereas spots, zonal 
bands, and sateUites all produce deviations of order ~ l%5p. These semi-analytic estimates 
are supported by more detailed numerical simulations. We also find that the light curve 
of a planet with the surface brightness profile expected from Lambert scattering deviates 
from that of a uniform source by ~ 10%^p. This affords the possibility of probing the phys- 
ical processes of the atmospheres of distant extrasolar planets by constraining their surface 
brightness profiles, and therefore the scattering properties of their constituent particles. 

We assessed the detectability of spots, rings, satellites, and bands with current and 
future telescopes. We found that, for reasonable assumptions and lOm-class telescopes, the 
planetary deviation will have a signal-to-noise of ~ 15, a ring system will only be marginally 
detectable with a signal-to-noise of ~ 6, and all other features will be completely unde- 
tectable. For 30m-class or larger telescopes, rings should be easily detectable. The detection 
of the non-uniform nature of the planetary surface brightness profile arising from Lambert 
scattering requires lOOm-class telescopes for bare detection. Spots, satellites and zonal bands 
are essentially undetectable for even the largest telescopes apertures. 
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Fig. 1. — An example of the caustic structure (thick sohd curve) produced by a binary lens. 
In this case, the system has equal mass components separated by one Einstein ring radius. 
The coordinates are centered at the midpoint of the binary, and all lengths are normalized by 
the Einstein ring radius. The solid and dashed straight lines are the two source trajectories 
considered in the numerical simulations. The insets are the details of the regions around the 
second caustic crossings of these trajectories. 
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Fig. 2. — Top: Caustic-crossing microlensing light curves from a planet with and without a 
satellite with radius equal to 20% of the radius of the planet, and albedo equal to 10% of 
the planet's albedo. The solid line shows the light curve from the planet only. The other 
lines show the magnification including a satellite with various relative positions. The events 
arise from the second caustic crossing of the dashed source trajectory through the binary- 
lens systems depicted in Fig. 1. Bottom: The additional magnification A^ga caused by the 
satellite alone. 
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Fig. 3. — Top: Caustic-crossing microlensing light curves from a planet with and without a 
face-on ring with different gaps between the planet's disk and the ring. The solid line shows 
the light curve from the planet only. The icons show the ring geometries for the corresponding 
hght curves with rings. The inner ring radius is R^i — 1.6 (dotted), 2.4 (short-dashed), and 
4.0i?p (long-dashed), respectively, where Rp is the radius of the planet. The outer ring radius 
is adjusted so that the projected area of the ring is equal to 7^{Rlut ~ ^in) — 10i?p. The 
ring has an albedo equal to 15% of the albedo of the planet. The lens system and the source 
trajectory responsible for the events are the same as for the light curves presented in Fig. 2. 
Bottom: The additional magnification AAj-i caused by the ring alone. 
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Fig. 4. — Top: Caustic-crossing microlensing light curves from a planet with inclined rings of 
different widths. The rings have a common inner radius of i?j„ = 2.0 Rp, but different outer 
ring radii of Rout = 2.6, 3.0, and 4.0 Rp. The ring has an inclination of i = 75°, and albedo 
relative to the planet of 15%. The lens system and the source trajectory responsible for 
the events are the same as for the light curves presented in Fig. 2. Bottom: The additional 
magnification AAy.^ caused by the ring alone. 
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Fig. 5. — Top: Caustic-crossing microlcnsing light curves from a planet with inclined rings 
with different albedos. The total magnification is shown for a planet with rings of equal size, 
[inner and outer ring radii of {Rin,Rout) = (2.0, 2.6)i?p], but relative albedos of = 0.075 
(solid), 0.3 (dotted), and 0.6 (dashed). The inclination of the rings is i — 75°. The lens 
system and the source trajectory responsible for the events are the same as for the light 
curves presented in Fig. 2. Bottom: The additional magnification AA^i caused by the ring 
alone. 
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Fig. 6. — Top: Caustic-crossing microlensing light curves from a planet with a ring, where 
the planet's shadow is cast on the ring. Note that, due to the specific geometry of the planet, 
host star, and observer, the planet is at quarter phase. Each ring has inner and outer ring 
radii of {Rin, Rout) = (2.4, 3.0)i?p, is seen face-on, and has an albedo relative to the planet 
of 15%. The lens system and the source trajectory responsible for the events arc the same 
as for the light curves presented in Fig. 2. For all line types, the heavier curve is for the 
planet only, while the lighter curve is for the planet and the ring. Bottom: The additional 
magnification A Ah caused by the ring alone. 
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Fig. 7. — Top: Caustic-crossing microlensing light curves from a planet with and without a 
spot with radius equal to 20% of the radius of the planet, and albedo equal to 80% of the 
planet's albedo. The sohd hne shows the hght curve from the planet only, the magnitude of 
the deviation due to the planet. The other lines show the magnifications including a spot 
with various relative positions. The lens system and the source trajectory responsible for 
the events are the same as for the hght curves presented in Fig. 2. Bottom: The deviation 
A^sp due to the spot. 
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Fig. 8. — Top: Caustic-crossing microlensing light curves from a planet with and without 
surface albedo variations caused by zonal bands. The zonal bands are parallel with the 
planet's equator and cause alternate variations of the surface brightness, and the planet is 
seen pole-on. The albedos of the bright and dark regions differ by 20%, and the mean surface 
brightness of each model has been adjusted in order to match that from the uniform planet. 
The events arise from the second caustic crossing of the solid source trajectory through the 
binary- lens systems depicted in Fig. 1. Bottom: The deviation AAba caused by the zonal 
bands. 
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Fig. 9. — Top: Caustic-crossing microlensing light curves from a planet with and without 
zonal bands. The zonal bands are parallel with the planet's equator and cause alternate 
variations of the surface brightness, and the inclination of the planet is 90°. The albedos 
of the bright and dark regions differ by 20%. The solid curve is the lensing hght curve 
corresponding to a planet with uniform surface brightness. The lens system and the source 
trajectory responsible for the events are the same as for the light curves presented in Fig. 8. 
Bottom: The deviation A^lba caused by the zonal bands. 
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Fig. 10. — Top: Caustic-crossing microlensing lensing light curves from a planet with the 
surface brightness profile expected from Lambert scattering (dotted) , and a uniform surface 
brightness profile (solid). The lens system and the source trajectory responsible for the 
events are the same as for the fight curves presented in Fig. 8. Bottom: The deviation A^La 
in the magnification of a Lambert sphere from a source with uniform surface brightness. 



